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ABSTRACT 



Oj' We study the possible phenomenology of a three-family Pati-Salam model constructed from 

intersecting D6-branes in Type II A string theory on the T®/(Z2 x Z2) orientifold with 
some desirable semi-realistic features. In the model, tree-level gauge coupling unification 
is achieved automatically at the string scale, and the gauge symmetry may be broken to 
the Standard Model (SM) close to the string scale. The small number of extra chiral exotic 
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C^ ' states in the model may be decoupled via the Higgs mechanism and strong dynamics. We 

o : 

calculate the possible supersymmetry breaking soft terms and the corresponding low-energy 

supersymmetric particle spectra which may potentially be tested at the Large Hadron Col- 

f^ . lider (LHC). We find that for the viable regions of the parameter space the lightest CP-even 

Higgs boson mass usually satisfies mn < 120 GeV, and the observed dark matter density 

^S ' may be generated. Finally, we find that it is possible to obtain correct SM quark masses 

5t 1 s-nd mixings, and the tau lepton mass at the unification scale. Additionally, neutrino masses 

and mixings may be generated via the seesaw mechanism. Mechanisms to stabilize the open 

and closed-string moduli, which are necessary for the model to be truly viable and to make 

definite predictions are discussed. 



I. INTRODUCTION 

Although string theory has long teased us with her power to encompass all known physical 
phenomena in a complete mathematical structure, an actual worked out example is still lacking. 
Indeed, the major problem of string phenomenology is to construct at least one realistic model 
with all moduli stabilized, which completely describes known particle physics as well as potentially 
being predictive of unknown phenomena. With the dawn of the Large Hadron Collider (LHC) era, 
new discoveries will hopefully be upon us. In particular, supersymmetry is expected to be found 
as well as the Higgs states required to break the electroweak symmetry. Therefore, it is highly 
desirable to have complete, concrete models derived from string theory which are able to make 
predictions for the superpartner spectra, as well as describing currently known particle physics. 

In the old days of string phenomenology, model builders were primarily focused on weakly 
coupled heterotic string theory. However, with the advent of the second string revolution, D- 
branes [l| have created new interest in Type I and II compactifications. In particular. Type IIA 
orientifolds with intersecting D6-branes, where the chiral fermions arise at the intersections of 
D6-branes in the internal space [4], with T-dual Type IIB description in terms of magnetized 
D-branes [3j, have shown great promise during the last few years. Indeed, intersecting D-brane 
configurations provide promising setups which may accommodate semi-realistic features of low- 
energy physics. Given this, it is an interesting question to see how far one can get from a particular 
string compactification to reproducing the finer details of the Standard Model as a low-energy 
effective field theory. 

In order to construct globally consistent vacua with intersecting D-branes, conditions must 
be imposed which strongly constrain the models. In particular, all Ramond-Ramond (RR) tad- 
poles must be cancelled and K-theory [4] conditions for cancelling the nontrivial Z2 anomally 
also must be imposed. Despite the clear benefits of supersymmetry, there have been many three- 
family standard-like models and Grand Unified Theories (GUT) constructed on Type IIA orien- 
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tifolds [a, la, lZ| which are not supersymmetric. Although these models are globally consistent, they 
are generally plagued by the gauge hierarchy problem and vacuum instability which arises from 
uncancelled Neveu-Schwarz-Neveu-Schwarz (NSNS) tadpoles. Later, semi-realistic supersymmetric 
Standard-like, Pati-Salam, unfli ppe d SU(5) as well as flipped SU(5 ) n iodels in Type IIA theory 



on T6/(Z2 X Z2) [S 



y, y, y, y, y, y, y, y and t6/(Z2 x z'^) y, y, y 



orientifolds were 
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eventually constructed, and some of their phenomenological consequences studied [l9|, l20(]. Other 
supersymmetric constructions in Type IIA theory on different orientifold backgrounds have also 



been discussed 



2l|. Nonperturbative D-instanton effects have also been receiving much attention 



of late, and may play an important role 22] [23|] [2j] 

In addition to satisfying the above consistency conditions, all open and closed-string moduli 
must be stabilized in order to obtain an actual vacuum. Unstabilized moduli are manifest in the 
low-energy theory as massless scalar fields, which are clearly in conflict with observations. Given a 
concrete string model, the low-energy observables such as particle couplings and resulting masses 
are functions of the open and closed string moduli. In a fully realistic model, these moduli must 
therefore be stabilized and given sufficiently large masses to meet the astrophysical/cosmological 
and collider physics constraints on additional scalar fields. Although satisfying the conditions for 
N = 1 supersymmetry in Type IIA (IIB) fixes the complex structure (Kahler) moduli in these 
models, the Kahler (complex structure) and open-string moduli generally remain unfixed. To 
stabilize some of these moduli, supergravity three-form fluxes 26l] and geometric flux es [27j| were 



introduced and flux models on Type II orientifolds have been constructed 
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31|, 



33, 



35l . |36|, |37l]. Models where the D-branes wrap rigid cycles, thus freezing the open-string moduli 



have also been studied 



161, 



13, 
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Despite substantial progress, there have been other roadblocks in constructing phenomeno- 
logically realistic intersecting D-brane models, besides the usual problem of moduli stabilization. 
Unlike heterotic models, the gauge couplings are not automatically unifled. Additionally, there 
has been a rank one problem in the Standard Model (SM) fermion Yukawa matrices, preventing 
the generation of mass for the first two generations of quarks and leptons. For the case of toroidal 
orientifold compactifications, this can be traced to the fact that not all of the Standard Model 
fermions are localized at intersections on the same torus. However, one example of an intersecting 
D6-brane model on Type IIA T^/(Z2 x Z2) orientifold has recently been discovered where these 
problems may be solved Q,Q. Thus, this particular model may be a step forward to obtaining 
realistic phenomenology from string theory. Indeed, as we recently discussed [38], it is possible 
within the moduli space of this model to obtain correct quark mass matrices and mixings, the tau 
lepton mass, and to generate naturally small neutrino masses via the seesaw mechanism. Further- 
more, it is possible to generically study the soft supersymmetry breaking terms, from which can 
be calculated the supersymmetric partner spectra, the Higgs masses, and the resulting neutralino 
relic density. 

This paper is organized as follows. First, we will briefly review the intersecting D6-brane model 
on Type IIA T^/(Z2 x Z2) orientifold which we are studying and discuss its basic features. We 
then discuss the low-energy effective action, and show that the tree-level gauge couplings are unifled 



near the string scale. We also find that the hidden sector gauge groups will become confining at a 
high energy scale, thus decoupling chiral exotics present in the model. Next, we study the possible 
low-energy superpartner spectra which may arise. We also calculate the Yukawa couplings for 
quarks and leptons in this model, and show that we may obtain the correct quark masses and 
mixings and the tau lepton mass for specific choices of the open and closed string-moduli VEVs. 
We should emphasize that for the present work, we will not focus on the moduli stabilization 
problem, as our goal is only to explore the possible phenemological characteristics of the model. 
However, we do comment on this issue and discuss how it may potentially be solved for this model. 
We also should note that models with an equivalent observable sector have been constructed in 
Type IIA and Type IIB theory as Ads and Minkowski flux vacua 



M, 



39[, so that the issue of 



closed-string moduli stablization has already been addressed to some extent. 

II. A D-BRANE MODEL WITH DESIRABLE SEMI-REALISTIC FEATURES 

In recent years, intersecting D-brane models have provided an exciting approach towards con- 
structing semi-realistic vacua. To summarize, D6 branes (in Type IIA) fill three-dimensional 
Minkowski space and wrap 3-cycles in the compactified manifold, with a stack of A^ branes having 
a gauge group U{N) (or U{N/2) in the case of T^ /{'L2 x Z2)) in its world volume. The 3-cycles 
wrapped by the D-branes will in general intersect multiple times in the internal space, resulting in 
a chiral fermion in the bifundamental representation localized at the intersection between different 
stacks. The multiplicity of such fermions is then given by the number of times the 3-cycles intersect. 
Due to orientifolding, for every stack of D6-branes we must also introduce its orientifold images. 
Thus, the D6-branes may also have intersections with the images of other stacks, also resulting in 
fermions in bifundamental representations. Each stack may also intersect its own images, resulting 
in chiral fermions in the symmetric and antisymmetric representations. The different types of rep- 
resentations that may be obtained for each type of intersection and their multiplicities are shown in 
Table HI In addition, there are constraints that must be satisfied for the consistency of the model, 
namely the requirement for Ramond-Ramond tadpole cancellation and to have a sprectrum with 
N = 1 supersymmetry. 

Intersecting D-brane configurations provide promising setups which may accommodate semi- 
realistic features of low-energy physics. Given this, it is an interesting question to see how far one 
can get from a particular string compactification to reproducing the finer details of the Standard 
Model as a low-energy effective field theory. There have been many consistent models studied. 



but only a small number have the proper structures to produce an acceptable phenomenology. 
A good candidate for a realistic model which may possess the proper structures was discussed 
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in Type IIA theory on the T^/fZo x Z2) orientifold. This background has been 



extensively studied and we refer the reader to 



9|] for reviews of the basic model building rules. 



We present the D6-brane configurations and intersection numbers of this model in Table IIIl and 
the resulting spectrum which is essentially that of a three-family Pati-Salam in Table IIIII 



We put the a', b, and c stacks of D6-branes on the top of each other on the third two torus, and 
as a result there are additional vector-like particles from N = 2 subsectors. 

TABLE I: General spectrum for intersecting D6-branes at generic angles, where laa' = ^2'^^'^' J|j_-^(ny^), 
and Ia06 — 2^^'''(~C'a^a + ll''^l''^l + ''^IPa'^^l + '^i'^a'a)' Morcover, M is the multiplicity, and as and a^ 
denote the symmetric and anti-symmetric representations of U{Na/2), respectively. 



Sector 


Representation 


aa 


U{Na/2) vector multiplet and 3 adjoint chiral multiplets 


ab + ba 


Mi^, f ) = /., = 2-'^ UliiKll - nlll) 


ah' + b'a 


M{^, f ) = lab' = -2-^ nliK^b + «X) 


aa -\- a a 


M{as) = \{Iaa' - ^laOe) ] M{aA) == ^{laa' + ^laOo) 



The anomalies from three global U{\)s of f^(4)c, U{2)l and C/(2)/j are cancelled by the Green- 
Schwarz mechanism, and the gauge fields of these U{l)s obtain masses via the linear Bf\F couplings. 
Thus, the effective gauge symmetry is SU{^)c x SU{2)l x 5C/(2)/j. In order to break the gauge 
symmetry, on the first torus, we split the a stack of D6-branes into ai and 02 stacks with 6 and 
2 D6-branes, respectively, and split the c stack of D6-branes into ci and C2 stacks with two D6- 
branes for each one, as shown in Figure [TJ In this way, the gauge symmetry is further broken to 
SU{^)c X SU{2)l X C/(l)/3^ X U{1)b-l- Moreover, the U{l)i^j^ x U{1)b~l gauge symmetry may 
be broken to U{1)y by giving vacuum expectation values (VEVs) to the vector-like particles with 
the quantum numbers (1, 1, 1/2, -1) and (1, 1, -1/2, Ij^under the SU{3)c x SU{2)l x C/(1)/3^ x 
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U{1)b~l gauge symmetry from a2c'i intersections 

Since the gauge couplings in the Minimal Supersymmetric Standard Model (MSSM) are unified 
at the GUT scale ~ 2.4 x 10^^ GeV, the additional exotic particles present in the model must 
necessarily become superheavy. To accomplish this it is first assumed that the USp{2)i and 



TABLE II: D6-brane configurations and intersection numbers for the model on Type IIA T^/Z2 x Z2 
orientifold. The complete gauge symmetry is [C/(4)c x U{2)l x C/(2)ij]obscrvabic x [C/S'p(2)'']hiddcn, the SM 
fermions and Higgs fields arise from the intersections on the first two-torus, and the complex structure 
parameters are 2xi — 6x2 — 3x3 — 6. 





f/(4)c x U{2)l X (7(2)« X t/5p(2)* 




iV 


(n\Zi) X {n^,f) X (n^P) 


Tig 


HA 
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b' 
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c' 
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a 
b 
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8 
4 
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(0,-l)x(l,l)x(l,l) 

(3,1) X (1,0) X (1,-1) 

(3,-l)x(0,l)x(l,-l) 



2 
-2 



-2 
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3 





-3 
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-1 


-1 
1 





3 



-3 




1 
2 
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2 
2 
2 
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(1,0) X (1,0) X (2,0) 
(1,0) X (0,-1) X (0,2) 
(0,-1) X (1,0) X (0,2) 
(0,-1) X (0,1) X (2,0) 


Xi = 3, X2 = 1, X3 = 2 
/3f = -3, PI = -3 
/3| = -3, /3f = -3 
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FIG. 1: Breaking of the effective gauge symmetry SU{A) x SU{2)l x SU{2)r down to S'C/(3)c x [/(2)l x 
U{^)i3R. X t/(l)_B-L via brane splitting. This process corresponds to giving a VEV to adjoint scalars, which 
arise as open-string moduli associated with the positions of stacks a and c in the internal space. 

USp{2)2 stacks of D6-branes lie on the top of each other on the first torus, so we have two pairs 
of vector-like particles XI2 with USp{2)i x USp{2)2 quantum numbers (2,2). These particles can 
break USp{2)i x USp{2)2 down to the diagonal USp{2)di2 near the string scale by obtaining VEVs, 
and then states arising from intersections al and a2 may obtain vector-like masses close to the 
string scale from superpotential terms of the form 



W C Xal Xa2 X\2. 
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TABLE III: The chiral and vector-like superficlds, and their quantum numbers under the gauge symmetry 
SU{A)c X SU{2)l X SU{2)ii x USp{2)i x USp{2)2 x USp{2)3 x USp{2)4,. 





Quantum Number 


Q4 


Q2L 


Q2R 


Field 


ah 
ac 
al 
o2 
62 
64 
cl 
c3 
6s 
6a 
cs 

CA 


3 X (4, 2, 1,1 
3 X (4,1,2,1 
1 X (4,1,1,2 
1 X (4, 1, 1, 1 
Ix (1,2, 1,1 
3 X (1,2,1,1 

1 X (1,1,2,2 
3x(l, 1,2,1 

2 X (1,3, 1,1 
2x(l,T,l,l 
2 X (1,1,3,1 
2 X (1,1,1,1 


1 
1 
1 
2 
2 




1 

-1 
1 

-1 










-1 




1 

-1 


2 
-2 






1 





-1 
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-2 
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Fl{Ql,Ll) 
Fr{Qr,Lr) 
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Sh 
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3 X (4,2,1,1 
3 X (4, 2, 1,1 






1 
-1 


1 
-1 








ad 


3 X (4,1,2,1 
3 X (4, 1,2,1 






1 
-1 





1 
-1 




he 


6 X (1,2,2,1 
6 X (1,2,2,1 










1 
-1 


-1 
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H'u, H'a 



where we neglect the couplings of order one. Moreover, we assume that the T^ and S'^ obtain 
VEVs near the string scale, and their VEVs satisfy the D-flatness of U{1)r. We also assume that 
there exist various suitable high-dimensional operators in the effective theory, and thus the adjoint 
chiral superfields may obtain GUT-scale masses via these operators. With Tj^ and 5*)^, we can give 
GUT-scale masses to the particles from the intersections cl, c3, and cs via the supepotential: 

W C 5)jXa X,i + TI,Xi, X^, + ^5),5^T|4. (2) 

The beta function for USp{2)di2 is —4 and the gauge coupling for USp{2)di2 will become strongly 
coupled around 5 x 10"*^^ GeV, and then we can give 5 x 10"*^^ GeV scale VEVs to S\^ and preserve 
the D-flatness of U{1)l. The remaining states may also obtain intermediate scale masses via the 
operators 



W C Xb2 X 1)2 SI H — X{,4 Xm- 

Mx 



(3) 



To have one pair of light Higgs doublets, it is necessary to fine-tune the mixing parameters of the 
Higgs doublets. In particular, the /x term and the right-handed neutrino masses may be generated 
via the following high-dimensional operators 



ijkl rnnkl 



(4) 



ijkl 



where y^ and y]^"^ are Yukawa couplings, and Mst is the string scale. Thus, the ^ term is TeV 
scale and the right-handed neutrino masses can be in the range 10^''"^^ GeV for y^fl ~ 1 and 

III. THE N = 1 LOW-ENERGY EFFECTIVE ACTION 



In building a concrete string model which may be testable, it is not enough to simply reproduce 
the matter and gauge symmetry of the known low-energy particle states in the Standard Model. 
It is also necessary to make predictions regarding the superpartner spectra and Higgs masses. If 
supersymmetry exists as expected and is softy broken, then it is possible to calculate the soft SUSY 
breaking terms, which determine the low energy sparticle spectra. Furthermore, if the neutralino 
is the lightest supersymmetric particle (LSP), then it is expected to make up a large fraction of 



id, 



4l|, and this is calculable 



the observed dark matter density, 0.0945 < ^h? < 0.1287 at 2a 

from the soft terms. Ideally, one would also like to be able to calculate the Yukawa couplings for 

the known quarks and leptons, and be able to reproduce their masses and mixings. 

To discuss the low-energy phenomenology we start from the low-energy effective action. From 
the effective scalar potential it is possible to study the stability 
19, 



M 



Yukawa coup 
investigated 



mg s 



47 



, the tree-level gauge couplings 



4J], gauge threshold corrections 45(, and gauge coupling unification 46|]. The effective 



48l |. matter field Kahler metric and soft-SUSY breaking terms have also been 



49( 1 . A more detailed discussion of the Kahler metric and string scattering of gauge. 



though 



28 



m, 



turningon Type IIB 3-form fluxes 



3d, 



31 



m 



33( 1 , there are additional 



matter, and moduli fields has been performed in [5d]. A 

can break supersymmetry from the closed string sector 

terms in the superpotential generated by the fluxes and there is currently no satisfactory model 

which incorporates this. Thus, we do not consider this option in the present work. In principle, 

it should be possible to specify the exact mechanism by which supersymmetry is broken, and thus 

to make very specific predictions. However, for the present work, we will adopt a parametrization 

of the SUSY breaking so that we can study it generically. 

The M = 1 supergravity action depends upon three functions, the holomorphic gauge kinetic 



function, /, Kahler potential K, and the superpotential W. Each of these will in turn depend 
upon the moduli fields which describe the background upon which the model is constructed. The 



holomorphic gauge kinetic function for a D6-brane wrapping a calibrated three-cyce is given by 34 1 



h 



2ttB 



n, 



Reie-'^^n 



'Up 



In terms of the three-cycle wrapped by the stack of branes, we have 

3 
^3 



n,. 



i=l 



from which it follows that 



(Up Tip Tip S 



Tl r> I r> t p 11 7T/p6ptr; 



Op hp ip 



plpip ' 



Ti p L p L p U ^ 



(5) 



(6) 



(7) 



^^^^..^..^..^^ 2(/32+/33) 2iP^+M 2('3i+/32) '' 

where Kp = 1 for SU{Np) and up = 2 for USp{2Np) or S0{2Np) gauge groups and where we 
use the s and u moduli in the supergravity basis. In the string theory basis, we have the dilaton 
S, three Kahler moduli T*, and three complex structure moduli C/* [50]. These are related to the 
corresponding moduli in the supergravity basis by 

Vlmt/iImC/2lmC/3\ 



Re(s) 

Re(u^) 
Re(t^) 




imuw^l 



ImW 



Im U'^ Im [/' 



jjkjji 



W 



{j,k,l) = (1,2,3) 



(8) 



and (/>4 is the four-dimensional dilaton. To second order in the string matter fields, the Kahler 
potential is given by 



K{M, M, C, C) = k{M, M)+ Y. ^c.c, (M, M)CiCj + 

untwisted i,j 

Y, Kc^Co(M,M)CeCe. 



(9) 



twisted 9 
The untwisted moduli Cj, Cj are light, non-chiral scalars from the field theory point of view, 
associated with the D-brane positions and Wilson lines. These fields are not observed in the MSSM, 
and if they were present in the low energy spectra may disrupt the gauge coupling unification. 
Clearly, these fields must get a large mass through some mechanism. One way to accomplish this 
is to require the D-branes to wrap rigid cycles, which freezes the open string moduli [171 ]. However, 
there are no rigid cycles available on T^ /{'L2 x ^2) without discrete torsion, thus we will assume 
that the open-string moduli become massive via high-dimensional operators. 
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For twisted moduli arising from strings stretching between stacks P and Q, we have ^ • 9pQ = 0, 
where 6pQ = (^q~ ^P ^^ ^^^ angle between the cycles wrapped by the stacks of branes P and Q on 
the j*^ torus respectively. Then, for the Kahler metric in Type IIA theory we find the following 
two cases: 



(^'po < 0, 0%o > 0, ^'po > 



'PQ 



'PQ 



PQ 



K 



PQ 



e^"e 



a,ie(2-y:U^pq) 



ne{ 



PQ' 



r(i - e^o) 



\ni+e'pQ)\ no%Q) \ ne'pQ 



Til 



'PQ> 



-k\-i+eir,{4 I ji\-i+e' 



{t^ + t^ypQ (r + t'^y+'^pQ (t' + 1 



I 
PQ. 



(10) 



o'po < 0, e'^o < 0, 9'po > 



'PQ 



'PQ 



PQ 



K 



PQ 



e^"e 



'4^7i5(2+E 



j=l °PQ 



ni+o%Q) 



ne\ 



PQI 



\ n-o'pQ) \ r(-e^Q) \\ r(i - 0^Q 



A , Tl\ 



{t^+t^y'-''pQ{r + i''y'~''pQ{f + i 



^PQ 



(11) 



For branes which are parallel on at least one torus, giving rise to non-chiral matter in bifunda- 
mental representations (for example, the Higgs doublets), the Kahler metric is 



K = {{s + s){t' + f){t^ + f){u^ + u^))-^/\ 



:i2) 



The superpotential is given by 

W = W + ]^^i^p{M)C^C^ + ^y(M)a/37C"'3T + • • • (13) 

while the minimum of the F part of the tree-level supergravity scalar potential V is given by 



N 1 



nM 



V{M, M) = e'^iGMK^'^'GN - 3) = (F'^RmmF'''' - 3e 



(14) 



where Gm = QmG and Knm = QnOmK, K is inverse of K^m, and the auxiliary fields F 
are given by 



M 



pM ^ ^G/2j^MLq^_ 



(15) 



Supersymmetry is broken when some of the F-terms of the hidden sector fields M acquire VEVs. 
This then results in soft terms being generated in the observable sector. For simplicity, it is assumed 
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in this analysis that the L'-term does not contribute (see 51[) to the SUSY breaking. Then the 



goldstino is included by the gravitino via the superHiggs effect. The gravitino then obtains a mass 

m3/2 = e^/2^ (16) 

which we will take to be ~ 1 TeV in the following. The normalized gaugino mass parameters, 
scalar mass-squared parameters, and trilinear parameters respectively may be given in terms of 
the Kahler potential, the gauge kinetic function, and the superpotential as 



m 



Iq = {ml/^ + Vo)-Y.F^F''dMdNlog{kpQ), 



M,N 



^PQR = F 



Km + dMlog{YpQR) - dMlog{KpQKQRKRp) 



where Km is the Kahler metric appropriate for branes which are parallel on at least one torus, i.e. 
involving non-chiral matter. 

The above formulas for the soft terms depend on the Yukawa couplings, via the superpotential. 
An important consideration is whether or not this should cause any modification to the low- 
energy spectrum. However, this turns out not to be the case since the Yukawas in the soft term 
formulas are not the same as the physical Yukawas, which arise from world-sheet instantons and 
are proportional to exp{—A), where A is the world-sheet area of the triangles formed by a triplet of 
intersections at which the Standard Model fields are localized. As we shall see in a later section, the 
physical Yukawa couplings in Type IIA depend on the Kahler moduli and the open-string moduli. 
This ensures that the Yukawa couplings present in the soft terms do not depend on either the 
complex-structure moduli or dilaton (in the supergravity basis). Thus, the Yukawa couplings will 
not affect the low-energy spectrum in the case of n-moduli dominant and mixed u and s dominant 
supersymmetry breaking. 

To determine the SUSY soft breaking parameters, and therefore the spectra of the models, we 
introduce the VEVs of the auxiliary fields Eq. (jlSp for the dilaton, complex and Kahler moduli 



m- 



F' = 2\/3Cm3/2Re(s)e,e-*^% 

F^'^'^y = 2V3Cmy2(.Mu')Q'^e-''^^ + Re(f )e*e-^^«*). (18) 

The factors 7^ and 7^ are the CP violating phases of the moduli, while the constant C is given by 



C^ = l + — ^. (19) 

^3/2 
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The goldstino is included in the gravitino by @s in S field space, and 6j parameterize the goldstino 
direction in [/* space, where X^(|0"P + |0iP) + |0sP = 1- The goldstino angle @s determines the 
degree to which SUSY breaking is being dominated by the dilaton s and/or complex structure (u*) 
and Kahler (t*) moduli. As suggested earlier, we will not consider the case of t-moduli dominant 
supersymmetry breaking as in this case, the soft terms are not independent of the Yukawa couplings. 

IV. GAUGE COUPLING UNIFICATION 

The MSSM predicts the unification of the three gauge couplings at an energy ~ 2.4 x 10^^ GeV. 
In intersecting D-brane models, the gauge groups arise from different stacks of branes, and so they 
will not generally have the same volume in the compactified space. Thus, the gauge couplings 
are not automatically unified, in contrast to heterotic models. For branes wrapping cycles not 
invariant under i}R, the holomorphic gauge kinetic function for a D6 brane stack P is given by 
Eq. d?!). where n* and s are the complex structure moduli and dilaton in the supergravity basis. 

The gauge coupling constant associated with a stack P is given by 

9m, = |Re(/p)|. (20) 

Thus, for the model under study the SU{3) holomorphic gauge function is identified with stack al 
and the SU{2) holomorphic gauge function with stack b. The Qy holomorphic gauge function is 
then given by taking a linear combination of the holomorphic gauge functions from all the stacks. 
Note that we have absorbed a factor of 1/2 in the definition of Qy so that the electric charge is 
given by Qem = ^3 + Qy- In this way, it is found [53] that 

Iy = -^fal + T,h2 + Tt/cI + o/c2- (21) 

5 2 2 2 

Recalling that the complex structure moduli C/* are obtained from the supersymmetry conditions, 

we have for the present model 

[/I = 3i , t/2 = 2 , [/3 = -1 + i . (22) 

Thus, we find that the tree-level MSSM gauge couplings will be automatically unified at the 
string scale 



2 2 5 2 

9s=9w = ^9y 



V6 ' 



(23) 



2tt 4 

Even though the gauge couplings are unified, this does not fix the actual value of the couplings as 
these still depend upon the value taken by the four-dimensional dilaton cp^. In order for the gauge 
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couplings to have the value observed for the MSSM (^unification ^ 0.511), we must choose (j)^ = —3 
such that e^'^"' ~ 20, which fixes the string scale as 

Mst = T^^/^e'*"'Mpi w 2.1 X lO^'^ GeV, (24) 

where Mpi is the reduced Planck scale. 

It should be kept in mind that values given for the gauge couplings at the string scale are only 
the tree-level results. There are one-loop threshold corrections arising from the A^ = 1 and N = 2 
open string sectors [45] which may alter these results. In addition, there is exotic matter charged 
under both observable and hidden sector gauge groups, which are expected to pick up large masses, 
but could still affect the running of the gauge couplings. 

V. CONFINEMENT OF THE HIDDEN SECTOR FIELDS 

In addition to the matter content of the MSSM, there is also matter charged under the hidden 
sector USp{2) gauge groups. These states will generally have fractional electric charges, similar to 



the so-called 'cryptons' 



54|, 



55l . l56l . l57l ] . Obviously, no such matter is observed in the low-energy 
spectrum so these exotic states must receive a large mass. Such a mass may arise if the hidden 
sector gauge couplings are asymptotically free and become confining at some high ener gy. Indeed, 



in the present case we find that the /3-functions for the USp{2) groups are all negative 



1^, 



PuSp{2)i — PuSp{2)2 — PuSp(2)3 — PuSp{2)4, — ~3, (25) 

where we consider all of the chiral exotic particles present even though it is expected that these 
states will decouple as discussed previously. From the holomorphic gauge kinetic function, the 
gauge couplings are found to take the values 

^2 „2 



9uSp(2)3 ~ 9uSp{2)4 



1. 



at the string scale. We may then straightforwardly run these couplings to low-energy energy via 
the one- loop RGE equations, 

where we find that the couplings for the USp{2)i and USp{2)2 hidden sector groups will become 
strong at a scale ~ 3 • 10-*^^ GeV, while the couplings for the USp(2)^ and USp{2)i groups will 
become strong around ~ 7 • 10^ GeV as shown in Figure [2l 
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USp(2)„ 




n (GeV) 



H (GeV) 



FIG. 2: RGB running of the gauge coupling for USp{2)i/USp{2)2, USp{2)3/USp{2)i hidden sector gauge 
groups, which become confining at « 3 • 10^^ GeV and « 6 • lO'^ respectively. 

We should note that it is also possible to decouple the chiral exotic states in the manner discussed 
in section II. 

VI. SOFT TERMS AND SUPERPARTNER SPECTRA 

Next, we turn to our attention to the soft supersymmetry breaking terms at the GUT scale 



defined in Eq. 
simplicity, as in 



8|. In the present analysis, not all the F-terms of the moduli get VEVs for 
5^. As discussed earlier, we will assume that Fj^ = so that the soft terms 



m, 



have no dependence on the physical Yukawa couplings. Thus, we consider two cases: 

1. The u-moduli dominated SUSY breaking where both the cosmological constant Vq and the 
goldstino angle are set to zero, such that F^ = F^ = 0. 

2. The u and s-moduli SUSY breaking where the cosmological constant Vq = and F^ ^ 0. 

A. SUSY breaking with u-moduli dominance 

For this case we take 0^ = so that the F-terms are parameterized by the expression 



F"" = V3m3/2(u* + u')eie-''^\ 



(28) 
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where i = 1,2,3 and with ^ |0jp = 1- With this parametrization, the gaugino mass terms for a 
stack P may be written as 



Mp = J^'^ Y, (^^^''' ®i ^~'^' n^prn'^pm^pj {j, k, I) = (1, 2, 3). (29) 

The Bino mass parameter is a hnear combination of the gaugino mass for each stack, 



\ k=l k=l 



MY = ^y^cpMp (30) 

JY ^ 

where the the coefficients cp correspond to the hnear combination of U{1) factors which define the 

hypercharge, U{1)y = ^cpt/(l)p. 

For the trihnear parameters, we have 

3 
ApQR = -\/3m3/2^ 
i=i 

+^"^3/2036-*'^^ (31) 

where P,Q, and R label the stacks of branes whose mutual intersections define the fields present 
in the corresponding trilinear coupling and the angle differences are defined as 

^PQ = (^Q - (^P- (32) 

We must be careful when dealing with cases where the angle difference is negative. Note for the 
present model, there is always either one or two of the 9pq which are negative. Let us define the 
parameter 

riPQ = sgn(Y[e'pQ), (33) 

i 

such that rjpQ = — 1 indicates that only one of the angle differences are negative while tjpq = +1 
indicates that two of the angle differences are negative. 

Finally, the squark and slepton (1/4 BPS) scalar mass-squared parameters are given by 



2 _ 2 
ITT-PQ — "^3/2 



1-3 y: 0™0n^"^^"'""""M%+E(^pr*(^PQ)+^pQ^pS*'(^pQ)) 

m,n=l \ jr=l 



.(34) 



The functions "^{Opq) = — -qq ^^ in the above formulas defined for tjpq = — 1 are 



if^PQ<0 : (35) 



Ide^.Q ^^' 2d9^,Q 



if^PQ>0 : 



^{e'po) = -^^ + \^ i^r(i - 9^pq) -\^ inr(^^pQ) + in(t^- + P), 
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and for ijpq = +1 are 



HO'pq) =^E + l^lnni + ey - i^ lnr(-0^po) - ln(t^- + P] 



iiepQ<0 : 

%) = 

iiOpQ>0 : 

'PQ 



^d9^PQ 



2d0^,,-— ^^«^ 



no'po) = iE + \-^ \^no'pQ) - \-^ inr(i - ey - in(t^- + P). 

atipf^ WpQ 



The function ^'(9pq) is just the derivative 



^'{9i 



d<i/{e{ 



PQl 



""'" <« 



ni,fc 



r.j,kl 



and 9pQ and 9pQ are defined [53] as 



(36) 



(37) 



^PQ 






r 1 ■ n -P 

[- 4^7 sin(27r6'-' ) J ^ when j = A; 

[^sm{27T9n]Q when j^k, 



e];g'=(n^ + n'=)(n' + n') 



52^i 



PQ 



du^du^ 



j^ [sm{A'K9^) + 4sm{2TT9^)]^ 



j^[sm{iTT9^)-4smi27r9^)] 

-j^[sm{A^9% 
lk[M^^9% 



Q 

p 
Q 



(38) 



when j = k = I 
when j ^ k = I 

when j = k^l or j = l^k 

when j ^ k ^ I ^ j. 

(39) 



Note that the only exphcit dependence of the soft terms on the u and s moduh is in the gaugino 
mass parameters. The trilinears and scalar mass-squared values depend explicitly only on the 
angles. However, there is an implicit dependence on the complex structure moduli via the angles 
made by each D-brane with respect to the orientifold planes. 

In contrast to heterotic string models, the gaugino and scalar masses are typically not universal 
in intersecting D-brane constructions, although in the present case, there is some partial universality 
of the scalar masses due to the Pati-Salam unification at the string scale. In particular, the trilinear 
A couplings are found to be equal to a universal parameter, ^o and the left-handed and right-handed 
squarks and sleptons respectively are degenerate. The Higgs states arise from the non-chiral sector 
due to the fact that stacks b, cl, and c2 are parallel on the third torus. The appropriate Kahler 
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metric for these states is given by Eq. (J12p . Thus, the Higgs scalar mass-squared values are found 
to be 

"i|/ = "^3/2(1-^103!'). (40) 

The complex structure moduli u^, and the four-dimensional dilaton (/>4 are fixed by the super- 
symmetry conditions and gauge coupling unification respectively. The Kahler modulus on the first 
torus t^ will be chosen to be consistent with the Yukawa mass matrices calculated in the next 
section. Thus, the free parameters which remain are 0i,2; sgn(03), f^, t^, the phases 7^, and 
the gravitino mass m^i2- In order to eliminate potential problems with electric dipole moments 
(EDM's) for the neutron and electron, we set 7j = 0. In addition, we set the Kahler moduli on 
the second and third tori equal to one another, Re(t^) = Re(t'^) = 0.5 and take the gravitino mass 
7713/2 ~ 1 TeV. Note that the soft terms only have a weak logarithmic dependence on the Kahler 
moduli. 

We constrain the parameter space such that neither the Higgs nor the squark and slepton scalar 
masses are tachyonic at the high scale, as well as imposing the unitary condition B^ + 02 + ©3 = 1. 
In particular, we require Q\ < 2/3, or equivalently 0f + ©^ > 1/3. 

To determine the soft terms and superpartner spectra at the low energy scale, we run the RGE's 
down from the high scale using the code SuSpect [60]. In principle, we should be able to determine 
tan /5, and the // and B parameters directly from the model. For the present construction, there 
is in fact a ^ parameter, whose real part corresponds geometrically to the separation between 
stacks b and c. In the absence of any effects which stabilize the open-string moduli, there are 
corresponding fiat directions in the moduli space. Indeed, the calculation of the Yukawa couplings 
in the next section will exploit this freedom. Thus, the effective //-term cannot be calculated until 
the moduli stablization issue has been addressed. Essentially the same considerations apply for a 
determination of tan /?, which depends upon the explicit values for the neutral component VEV of 
the Higgs fields, up to an over- all constant. The overall issue of moduli stablization is discussed in 
a later section. 

For the present work, we will fix these values via the requirement for electroweak symmetry 
breaking (EWSB), in a similar fashion to 



53]. We also choose /x > as favored by {g — 2) and 
take tan /? as a free parameter. We use the value for the top quark mass mt = 170.9 GeV. Then, 
knowing the low energy spectra, we can then determine the corresponding neutralino relic density. 
To calculate this, we use the code MicrOMEGAS [61] . Some points in the parameter space which 
give the observed dark matter density are shown in Table [IV] where we have fixed ©^ = 0. 
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TABLE IV: Supersymmetry breaking soft terms at the string scale and resulting neutralino relic density for 
some specific choices of goldstino angles, with G^ = 0. 



ei 


92 


63 


Mg (GeV) 


M^- (GeV) 


M^ (GeV) 


mu (GeV) 


rriL (GeV) 


niR (GeV) 


At (GeV) 


LSP 


nh^ 


-0.610 


0.290 


0.737 


889 


251 


422 


429 


963 


466 


676 


B 


0.115 


-0.610 


0.380 


0.695 


931 


329 


416 


524 


989 


478 


610 


B 


0.108 


-0.600 


0.470 


0.647 


967 


407 


411 


609 


1005 


492 


531 


B 


0.113 


-0.100 


0.870 


0.482 


1171 


753 


667 


806 


684 


521 


-312 


B 


0.118 


0.0400 


0.600 


0.799 


1211 


519 


920 


205 


510 


590 


—222 


B 


0.105 


0.110 


-0.570 


0.814 


211 


-493 


564 


74 


561 


631 


431 


B/H 


0.107 


0.150 


0.660 


0.736 


1209 


571 


944 


432 


529 


512 


-397 


B 


0.108 


0.010 


0.810 


0.586 


1209 


701 


793 


695 


603 


510 


-371 


B 


0.094 



TABLE V: Low energy supersymmetric particles and their masses (in GeV) for Qi 
63 = 0.737, and 6^ = 0. with tan /? = 46. 



-0.610, 62 = 0.290, 



h° 


//o 


A° 


H± 


9 


xf 


X^ 


X? 


xl 


116.89 


826.80 


826.81 


831.05 


1985 


192.3 


1115. 


173.3 


192.3 


xi 


xi 


ii 


t2 


Ul/Ci 


U2/C2 


&i 


&2 




-1113. 


1114 


1477 


1789 


1949 


1761 


1638 


1791 




di/Si 


di/h 


n 


r2 


J>T 


ei/fli 


62 //i2 


Ve/V^ 




1950 


1760 


189.0 


928.2 


919.5 


973.3 


488.0 


970.1 





A contour plot of the dark matter density is displayed in Fig [3] with 7713/2 = 1 TeV, tan /3 = 46 
and for 03 > 0. It can be seen that only small regions of the allowed parameter space can produce 
the observed dark matter density, which are indicated on the plot as dark bands. Regions which 
do not satisfy the constraints 



e? + e| + e^ = 1 



and 



1 



ef + e^,^ 



(41) 



(42) 



are indicated on the plots as the light gray shaded areas. The dark gray areas indicate regions of 
the parameter space where either the neutralino is not the LSP, LEP superpartner mass limits are 
not satisfied, or for which there is no RGE solution. The viable parameter space may be further 
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TABLE VI: Low energy supersymmetric particles and their masses (in GeV) for 81 = 0.04, Q2 — 0.60, 
63 = 0.799, and Qs = with tan f3 = 46. 



/J" 


H" 


r 


H^ 


ff 


X? 


xr 


X'l' 


X§ 


118.82 


1057.1 


1057.1 


1060.4 


2616. 


414.6 


1478. 


390.4 


414.6 


xi 


xi 


fi 


fa 


Ul/Ci 


U2/C2 


61 


62 




-1475. 


1477 


1955 


2114 


2318 


2335 


2071 


2184 




di/si 


d2/s2 


Tl 


r2 


!>T 


ei/Ai 


e2//i2 


J^e/t-M 




2319 


2330 


474.4 


689.6 


574.9 


620.1 


679.0 


615.1 





Relic Density 
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nh 



' 1.000 



-I 2.000 

1 3.000 

' 3.900 



FIG. 3: Contour map of the dark matter density for tan /3 = 46, 6^ > 0, Qg — and /i > 0. The areas in 
black denote regions where 0.09 < fth^ < 0.12 with a gravitino mass TO3/2 = 1 TeV. The light gray regions 
are excluded because they do not satisfy the constraints on the soft terms at high scale. The dark grey 
regions denote regions where either the neutralino is not the LSP, mass limits are not satisfied or for which 
there is no RGE solution. Regions inside the white contour satisfy the LEP limit, nih > 114 GeV. 



constrained by imposing the LEP limit on the lightest CP-even Higgs mass, mh > 114 GeV. 
Regions satisfying this bound are contained within the white contom' on the plot. Essentially, this 
rules out regions of the parameter space with a mixed Bino/Higgsino LSP as the Higgs mass for 
these regions of the parameter space is always below 114 GeV. A sampling of some of the SUSY 



spectra for some of the cases of Table IIVI are shown in Tables |V] and IVI[ 
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LSP - Stau Mass Difference 



LSP - Chargino Mass Difference 




I 





-0.8 -0.6 -0.4 -0.2 OCO 0.2 0.4 0.6 0.8 



FIG. 4: Contour map of the neiitralino-stau mass difference and neutrahno-chargino mass difference for 
tan (3 = 46, 6*3 > 0, 0s = and fj, > 0. The areas in black denote regions where 0.095 < flh^ < 0.125 
with a gravitino mass 777,3/2 — I TeV. The light gray regions are excluded because they do not satisfy the 
constraints on the soft terms at high scale. The dark grey regions denote regions where either the neutralino 
is not the LSP, mass limits are not satisfied or for which there is no RGE solution. The dark matter density 
within the allowed range coincides with the Xi "^ £^rid x1 x^/Xi X2 coanhiliation regions where the stau 
and/or chargino/next-to- lightest neutralino mass is slightly bigger than the lightest neutralino mass. 



It can be seen from Fig. [3]that regions of the corresponding parameter space correspond to cases 
where the lightest neutralino mass is very close to the either or both the lightest chargino/next- 
to-lightest neutralino Xi /X2 (which are mass degenerate) and light stau fi mass. In other words, 
the observed dark matter density is obtained close to the Xi '^i and/or x? Xi /Xi X2 coanhilation 
regimes. 

The lower bound on the gravitino mass may be estimated for fixed tan (3 by lowering 7^3/2 until 
there are no regions of the parameter space for which the Higgs mass satisfies the LEP limit. In 
this way, it is found that the lower bound satisfies 425 GeV < 7713/2 < 450 GeV for tan (3 = 46. 
This can be seen in Figs. [Sj where the region of the parameter space with a Higgs mass above 
114 GeV begins to shrink dramatically. For 771,3/2 = 450 GeV, only very small regions of the viable 
parameter space results in a Higgs mass above 114 GeV. For 7113/2 = 425 GeV, there is no region 
of the parameter space above this limit. Similar results hold for other values of tan (3. 
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Relic Density, m^^^ = 475 GeV, tan p = 46 



Relic density, m ^ 450 GeV, tan p ^46 
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FIG. 5: Contour maps of the relic density for tan /? — 46, ^3 > 0, B^ = and ^ > with a gravitino mass 
777,3^2 — 475 GeV and 7713/2 — 450 GeV respectively. The areas in black denote regions where 0.09 < fth^ < 
0.12. The light gray regions are excluded due to tachyonic scalar masses at the high scale. The dark grey 
regions denote regions where LEP superpartner mass limits are not satisfied or there is no RGE solution. 
Regions inside the white contour satisfy the LEP limit, rrih > 114 GeV. The viable region of the parameter 
space is much smaller for 777,3/2 — 450 GeV and dissappears for 7773/2 — 425 GeV, which constrains the 
minimum gravitino mass to be in the range 425 GeV < 7773/2 < 450 GeV. 

B. SUSY breaking via ?/-moduli and dilaton s 



Next, we allow the dilaton s to obtain a non-zero VEV as well as the ^-moduli. To do this, we 
parameterize the i^-terms as 



V3ms/2[{s + s)ese-'^^ + {u' + w^)e^e- 



-t^i 



(43) 



Clearly, this is a more complicated situation with a much larger set of values over which to scan. 
The formula for the gaugino mass associated with each stack now becomes 



Mp 



-VS: 



m 



3/2 



Re/^ 



Y^ Re(n^) Qj e"*^^ n^pm%m}p2-^^*'+^'^ 



+ esRe(s)e-*^On),n|,n^ 



(44) 



(j,/c,0 = (1,2,3). 



As before, the Bino mass parameter is a linear combination of the gaugino mass for each stack, and 
the coefficients corresponding to the linear combination of U{1) factors define the hypercharge. 
The trilinear parameters generalize as 



ApQR = -\/3m3/2^ 
V3 



e,e-^^ (i + {JZi%^{e%Q) -\) + {Jz&{oip) - \. 



k=l 



fe=l 



+^"i3/2(©3e-*^^+e,e- 



-«7s 



(45) 
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where 0o corresponds to G^ and there is a contribution from the dilaton via the Higgs (1/2 BPS) 
Kahler metric, which also gives an additional contribution to the Higgs scalar mass-squared values: 



2 _ 2 
rrijj^ — 771,3/2 



i-^(|e3|' + |e,|2) 



(46) 



Finally, the squark and slepton (1/4 BPS) scalar mass-squared parameters are given as before 



by 



2 2 

mpQ = m. 



■3/2 



1-3 j; e^e.e-(^™--")(%^ + ^(e];7^(e^PQ)+e^-erQ*'(0^Q)) 



_ _ , ,(47) 

m,n=0 \ jf=l 

where we now also include the Qg = @o in the sum. The functions ^{9pq) and ^'{Opq) are given 
as before by Eq. (f36l) and Eq. (f37l) . The terms associated with the complex moduli in ^pQ and 
^pg are the same as those in Eq. (j38p and Eq. (j39p . and the terms associated with the dilaton are 
given by 






Att 



sin(27re^), 



Q 



(48) 



e^^^^{u^ + u^){s + s) 



dH{ 



PQ 



du^ds 
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i- [sin AttO^} „ when j = k 



(49) 



'ih^ ['^™47r6'^]Q when j y^ k 



and 



e^Q^ ^{s + s){s + s) -^ = — [sinATre^ + 4sin(2vre^)] J , 



(50) 



where k,l ^ s. As before, the Gj parameters are constrained as Ylii=i ©? + ©s = 1- I'^ order to 
simplify the analysis, we fix Gs while varying G3, Gi, and Q2 subject to the unitary condition. 
Since there is now another free parameter, the possible parameter space is much larger than in the 
previous case. Some points in the parameter space which give the observed dark matter density 
are shown in Table IVlIl where we have fixed G^ = 0.40. 

We exhibit the relic density as a function of Gi and Q2 for the particular case with Q^ = 0.40 
and tan /? = 46 in Fig 6. 

In the figure, only the regions within the white contour satisfy the LEP limit on the Higgs mass. 
Essentially, this rules out regions of the parameter space with a mixed Bino/Higgsino LSP as the 
Higgs mass for these regions of the parameter space is always below 114 GeV. A sampling of some 
of the superpartner spectra for some of the cases of Table IVTIl are shown in Tables rVIIIl IIXI and 1X1 
below. The viable parameter space is quite large as can be seen from Fig. 7. 
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TABLE VII: Supersynimetry breaking soft terms at the string scale and resulting neutralino relic density 
for some specific choices of goldstino angles, with 6s — 0.40. 



ei 


62 


93 


Ma (GeV) 


M^ (GeV) 


Ms (GeV) 


ruH (GeV) 


ruL (GeV) 


niR (GeV) 


At (GeV) 


LSP 


nh^ 


-0.490 


0.640 


0.592 


932 


207 


344 


688 


832 


659 


154 


B 


0.112 


-0.450 


-0.45 


0.771 


181 


-736 


181 


327 


470 


447 


842 


B/H 


0.118 


0.140 


0.870 


0.473 


971 


407 


592 


815 


470 


562 


-732 


B 


0.114 


0.380 


-0.45 


0.808 


218 


-736 


649 


142 


564 


671 


-68 


B/H 


0.105 


0.480 


-0.41 


0.776 


220 


-701 


682 


312 


692 


638 


-213 


B/H 


0.112 


0.600 


-0.13 


0.789 


476 


-458 


856 


255 


765 


520 


-516 


B 


0.112 



Relic density, m^^^ = 1 TeV, tan p = 46, ©^ = 0.4 
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FIG. 6: Contour maps of the relic density for tan (3 = 46, 6s = 0.40 7^ 0, 6^3 = 0.75, and /i > with a 
gravitino mass 7713/2 = 1 TeV. The light gray regions are excluded because they do not satisfy the constraints 
on the soft terms at high scale. Regions inside the white contour satisfy the LEP limit, nih > 114 GeV. 



As before, it can be seen that regions of the corresponding parameter space correspond to cases 
where the the lightest neutralino mass is very close to the either or both the lightest chargino/next- 
to-lightest neutralino and light stau mass. Thus, the observed dark matter density is obtained close 
to the Xi "^1 and/or Xi X /Xi X2 coanhilation regimes. Similar results hold for other values of tan 
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TABLE VIII: Low energy supersymmetric particles and their masses (in GeV) for 9i 
-0.130, 63 = 0.789, and 6^ = 0.40 with tan (3 = 46 and 7713/2 = 1 TeV. 



-0.600, 62 



/l° 


H° 


A" 


H± 


9 


X? 


X^ 


X? 


X§ 


114.28 


600.19 


600.2 


606.28 


1135 


390 


774.8 


-363.2 


-390.0 


X§ 


xl 


ti 


t2 


Mi/ci 


W2/C2 


fei 


&2 




769.5 


-772.8 


856.5 


1133. 


1253. 


1105. 


985.4 


1126. 




di/si 


d2/s2 


Tl 


T2 


Vt 


ei//ii 


e2//i2 


Ve/V^i 




1255. 


1091. 


470.4 


796.9 


787.6 


834.2 


607.9 


830.5 





TABLE IX: Low energy supersymmetric particles and their masses (in GeV) for 0i 
83 = 0.592, and G^ == 0.40 with tan /? = 46 and mg/a = 1 TeV. 



-0.490, 62 = 0.640, 



/i" 


H° 


A" 


Ji-± 


S 


X? 


xt 


X? 


X^ 


117.82 


885.60 


885.61 


889.63 


2079 


155.1 


1110. 


139.1 


155.1 


xl 


xS 


ii 


t2 


Ul/Ci 


U2/C2 


&1 


62 




-1107 


1109 


1573. 


1774. 


1961. 


1900. 


1731 


1793. 




di/Si 


d2/s2 


fi 


r2 


i^T 


ei//ii 


62/^2 


VelVti 




1962. 


1900 


486.1 


793.6 


775.9 


838.1 


669.4 


834.5 





VII. YUKAWA COUPLINGS 



In addition to the fact that the SM fermions are replicated into three distinct generations, the 
different generations exhibit an intricate pattern of mass hierarchies. At present, there has been no 
satisfactory explanation for this. However, in addition to the replication of chirality, intersecting 
D-brane models may naturally give rise to mass hierarchies and mixing. 

The Yukawa couplings in the intersecting D-brane worlds arise from open string world-sheet 
instantons that connect three D-brane intersections [62]. For a given triplet of intersections, the 
minimal action world-sheets which contribute to the Yukawa coupling are weighted by a factor 
exp{—Aabc), where Aabc is the world-sheet area of the triangles bounded by the branes a,b, and 
c. Since there are several possible triangles with different areas, mass hierarchies may arise. One 
may also see that the Yukawa couplings depend on both the D-brane positions in the internal 
space, as well as on the geometry of the underlying compact manifold. Effectively, these quantities 
are parameterized by open and closed-string moduli VEVs. Thus, in order to make any definitive 
predictions for the Yukawa couplings, these moduli must be stabilized. 



25 



TABLE X: Low energy supersymmetric particles and their masses (in GeV) for 81 = 0.140, 82 = 0.870, 
63 = 0.473, and 6^ = 0.40 with tan (3 = 46 and 7713/3 = 1 TeV. 



h° 


H" 


A" 


H^ 


ff 


xf 


xt 


X? 


X^ 


118.02 


911.27 


911.17 


915.16 


2136 


323.6 


1138. 


247.4 


323.6 


X§ 


xl 


fi 


t2 


Ul/Cl 


U2/C2 


&i 


&2 




-1134 


1137 


1515. 


1691. 


1911. 


1933. 


1639 


1762. 




di/si 


dijai 


T\ 


T2 


i>T 


ei//ii 


62 //i2 


VejVp. 




1912. 


1931 


267.6 


517.4 


425.9 


540.1 


602.7 


534.4 






FIG. 7: The viable parameter space with W? < 0.125 for 8^ > 0. 



Despite the promise of intersecting D-brane models in explaining the fermion mass hierarchies 
and mixings of the Standard Model, they have typically been plagued by a rank-1 problem in 
the Yukawa mass matrices. This can be traced to the fact that not all of the Standard Model 
fermions are localized at intersections on the same torus (in the case of toroidal orientifold com- 
pactifications) . The general flavour structure and selection rules for intersecting D-brane models 



as been investigated in [631 . 



64( 1 However, the intersections are all on the first torus for the model 



of the previous section. Thus the resulting Yukawa mass matrices do not have a rank-1 problem. 
In the following, we will explore the moduli space of this model in order to see if there are any 
solutions which may give rise to realistic Yukawa textures, following the analysis of [47], without 
first addressing the issue of moduli stabilization. Our goal for the present work is simply to see if 
the observed fermion mass hierarchies may arise in this model and identify the points in the moduli 



26 



space where this may happen. We discuss how the moduh might be fixed in a later section. 

A. The general form 

We start by considering D6 branes in Type IIA, where the D6-branes wrap 3-cycles which 
intersect at angles on a compact manifold T^ = T^ x T^ x T^. For simplicity, we consider three 
stacks of D-branes wrapping on a two-torus. The 3-cycIes wrapped by the D-branes can be written 
in a vector form in terms of the wrapping numbers: 

[Ua] = na[a] + ma[b] : Za = R{n + Uma) ■ Xa, 
[lib] = nb[a] + nibib] : Zb = R{nb + Unib) ■ Xb, 
[nj = nc[a] + mc[b] : Zc = R{nc + UrUc) ■ Xc, (51) 

where U is the complex structure parameter of the torus and x € R arbitrary numbers. The 
Yukawa coupling involving branes a, b, and c recieves a contribution from the areas of the triangles 
bounded by the triplet of D-branes. To ensure that the triplet of intersections actually form a 
triangle, we must impose the closer condition 



47] 



Za + Zb + Zc = 0. (52) 

The wrapping numbers are all integers due to the quantization conditions, so by translating Eq. 
()52p into the Diophantine equation, the solution is found to be 

_ he 

Xb = —rX, X = Xq + I, Xq G H, I £ Z, 
a 

Xc ^ ~r'^i v5"jj 

where lab is the intersection number, and d = g.c.d.{Iab,Ibc,Ica) is the greatest common divisor of 
the intersection numbers. The parameter / indexes the different points in the covering space C but 
the same points in the lattice of T"^ of the triangles. The quantity xq is dependent on the different 
intersection points of two of the D-branes which are indexed by 

i = 0,l,--- ,(l/a6l-l), 

j = 0,l,--- ,{\hc\-l), 

k = 0,l,--- ,i\Ica\-l), (54) 
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therefore xq can be written as 



xo{^:J,k) = — + -;— + 



(55) 



It is not necessary to require that all branes intersect at the origin. If the position of the stacks 
is shifted by an amount ea, a = a,b,c, in clockwise directions of stack a by a length in units of 
^/IIIIqII on each torus, then xq will be modified as 



.. . ,x ^ , i , ^ , d{Iabec + Ica^b + lab^a) 

X(}{l,J,k) = — + ^— + -;— + 



(56) 



'ab ^ca ^bc ^ab^bc^ca 

With this parameterizion of xq, we can now calculate the areas of the triangles, by the area formula 
of vectors 



1 



(57) 



A{Za,Zb) = -i/|Za|2|z;,|2 - {ReZaZhY 

-^ Aijk{l) = \{2TlfA\IabIbcha\{^ + j- + — + e + If, 
^ i-ab ^ca ^bc 

where e is the total shift effect in Eq. ([56j) . 

J-ab^bc^ca 

and A is the Kahler structure of the torus, which is generally the area. By adding a real phase 
(^abc = s^S^il abhcica) for the full instanton contribution, the corresponding Yukawa coupling for 



(58) 



the three states localized at the intersections indexed by {i,j, k) is 



m 



Yijk = hquaabc ^ exp(- 



27ra' 



(59) 



where hgu is due to quantum correction as discussed in [48,]. The summation can be expressed in 
terms of a modular theta function for convenient numerical calculation. The real version of this 
theta function can be written as 



1} 



w = E 



^-nt{5+lf ^2lm{5+l)4> 



(60) 



/ez 



so after comparing the parameters, we have 



i j k 

4> = 0, 

t ^ 7 -'afe-'fec-'ca • 

a' 



(61) 
(62) 
(63) 
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a. B-field and Wilson lines 

The theta function above is constrained by its real property. However, t can be complex 
and (/) can be any number as an overall phase which can be given both a theoretical and 
phenomenological interpretation. 

If we turn on a B-field in the compact space T^, the string will not only couple to the metric 
but also to this B-field. Then the Kahler structure may be written in a complex form 



J = B + iA, 



(64) 



where the parameter t is replaced by a complex parameter k 

J 



a' 



l-'afe-'fec-'cal 



(65) 



We can also include Wilson lines around the compact directions that the D-branes wrap in this 
construction [471]. To avoid breaking any gauge symmetry, these Wilson lines are chosen up to a 
phase. If we consider the Yukawa coupling formed by D-branes a, b, and c, wrapping a different 
one-cycle with Wilson lines in the phases exp(27ri9a), exp(2iTi6f,), and exp(27ri^c) respectively, then 
the total phase is a linear combination of each phase weighted by the relative longitude of each 
segment, determined by the intersection points: 



2mXa9a J^-rTiXf,0i, 2TTiXc0c _ „2mx{Itcda+Ica0b + Iab9c) 



(66) 



Thus, including these two effects, we obtain a general complex theta function as 



where 



b. Other modification 



^ 



(^) = ^e™('^+')% 



27ri(<5+Z)</) 



(67) 



iez 



i j k 



lab 



he 



hcGa + IcaOh + lahOc, 
J , 



a' 



^ah^bc^ca\ ■ 



(68) 
(69) 
(70) 



In most of the (semi-)realistic models orientifold planes must be introduced to cancel the 
RR-tadpoles. As a result, there will be additional fields from the brane images coupling to the 
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ordinary branes fields as well as themselves. For example, for the triangle formed by branes o, 
b', and c, the Yukawa coupling will then depend on the parameters lab', h'a &iid lea, and the 
corresponding indicies are i', j, and k', where the primed indexes are independent of the unprimed 
ones. 

The other issue is the non-coprime cases. The three intersection numbers are not necessarily 
coprime, so we have to make sure we do not over-count the repeated parts. The constant d is 
defined as the g.c.d. of the intersection numbers and is introduced in the brane shift parameters. 
We must then modify the other parameters as well: 

, hcOa + Icadb + labdc 



d 
'J l^ab^bc-^cal 



(71) 
(72) 



a' d2 

There is one more constraint which is necessary to ensure non-zero Yukawa couplings: the 
triangles must be bounded by D-branes. Thus the intersection indexes i, j, and k cannot be 
arbitrary. They are required to satisfy [471] 



i + j + k = mod d. 



(73) 



There is one degree of freedom which may ease this constraint, which is an additional parameter 
in 6: M 



i j k s 

J-ab J-ca J^bc 0' 



(74) 



where s is a linear combination of i, j, and /c, and it is just a shift of counting the triangles since 
we have required {i, j. A;} to satisfy ([73]) . 

Therefore finally, we can write down a complete form for the Yukawa couplings for D6-branes 
wrapping on a full compact space T^ x T^ x T^ as 



^{ijk} = hquCJabc | J_ ^ 



r=l 






Arh 



(75) 



where 



Sir) 



ci-)) 



IrGZ 



^ 1 



(76) 
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with r = 1, 2, 3 denoting the three two-tori. The input parameters are given by 



T(r) I r(^) /-('-) /-('-) I 
•^ l-' ab -'fee -'ca I 



^M ^ I "':;;,,;"' • (77) 



c. Theta function with characters 

It can be comphcated to calculate the numerical value of the theta function defined in 
Eq. (|67|) . Thus, for simplicity the B-field will not be introduced in the following analysis. Then if 
we define J' = —iJ = A and so k' = —in for convenience, the theta function turns out to be 







" 






^ 


6 


(^') = y^ e-TK'(5+O%27ri(<5+O0 






[<p\ 


iGZ 


redefine 

> 


^ 


6 

_4>_ 


{K)=e-'"'^\^^'^'t>^^{Ti{(t> + iK 



(78) 



where ^3 is the Jacobi theta function of the third kind. A plot of the theta function is shown in 
Fig. H with = 0. 

B. Semi- Realistic Yukawa Textures 

(i) Mass Matrices 

As mentioned previously, the intersecting D-brane model with matter content shown in Table 
mi has several desirable semi-realistic features, namely three-generations of chiral SM fermions 
with a minimum of exotic matter, tree-level gauge coupling unification, and the fact that the three 
intersections required to form the disk diagrams for the Yukawa couplings all occur on the first torus 
as can be seen from Figure [9l Thus, in our analysis we will focus on the just the first torus, since 
the contribution from the other two tori will just give an over-all constant. This constant, which is 
different for the up-type and down-type quark, and charged lepton mass matrices, is unimportant 
for the present analysis since we will only obtain the mass ratios rather than the absolute fermion 
masses. 

As described in the previous analysis, the Pati-Salam gauge symmetry is broken to the Standard 
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FIG. 8: Theta function as a function of the parameters S and k, with (f> — 0. The maximum values occur 
for (5 = while the minimum values occur when 6 = 1/2. 

Model by a process which involves brane-splitting, 



O ^ Ol + 02, 



Cl +C2. 



(79) 



so that the Standard Model quarks and leptons arise from 



Fl{Ql,Ll)^Ql + L 
Fr{Qr, Lr) ^Ur + Dr + Er + N. 



(80) 



The Yukawa couplings for the quarks and leptons are then given by the superpotential 



(81) 



where it should be kept in mind that there are six vector pairs of Higgs multiplets, each of which 
may receive a VEV. 

For the model under consideration, the intersection numbers on each torus are given by 



'afe 



3 /(^) 
-3 /^) 
-6 /') 



-1 /(') 

-1 /(') 



-1, 



1, I, 



(3) 



fee 



0. 
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FIG. 9: Brane configuration for the three two-tori. The SM ferniion mass hierarchies primarily result from 
the intersections on the first torus. 



These intersection numbers are not coprinie, so that we need to have S^' = g.c.d.{I^^ , ^J , lea ) = 
3, d^'^' = 1. Note that we do not need d'^^ since the intersections do not form any triangles on the 
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third torus. Thus the parameters of the theta functions defined in Eq. [77| are given by 

g) _ ^^'^ j(^) kW eP - ej'^ - 2e^^ .(^) 
3 3 6 6 3 ' 

Si^) = -ef) + ef\ (82) 

(^i^) = -ef) + ef\ (83) 



.<.. = 5^ 



a' 

k(3) = 0, (84) 

where i = ... 2, j = ... 2, and k = ... 5, indexing the left-handed fermions, right-handed 
fermions, and Higgs fields respectively. For mathematical convenience we redefine the shift on each 
torus as 

6 

e^^)^-ef^ + ef\ (85) 

The intersection numbers on the second and third tori are either one or zero, so their contribution 
to the Yukawa couplings will be an over-all constant. For a given triplet of intersections to be 
connected by an instanton, the selection rule 

^(1) + j{i) + ^^(1) = mod 3 (86) 

should be satisfied. Then, the Yukawa coupling matrices will take the following form: 
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/ oooo 
ai2o 

\ 0210 J 









/ 0003 

ai23 

\ 0213 J 



, Y, 






O012 



\ 



O102 

\^ 0222/ 



/o 



0015 



\ 



^ / ao2i\ / 

oiii 

\0201 y 

^ / O024\ 

k=A- 0114 , n=5~ 0105 . (87) 

\0204 y y 0225/ 

By choosing a different linear function for s'^', some independent modes with non-zero eigen- 
values are available, which are listed below. 
(i) s(i) = 



Oooo — O102 — O204 — '& 



0210 — O012 — O114 



{} 



O12O = O222 = O024 = ^ 
O02I = O123 = O225 = ^ 
O2OI = O003 = Ol05 = '& 
Olll = O213 = O015 = ■& 



e(i) 

6(1) + 1 
0(1) 

^ 3 

0(1) 

e(i) + i 

0(1) 

e(i) + i 

" ^2 

0(1) 

" 6 

0(1) 



6J(1) 
(^)^A 



a' 



.6J(1), 



a' 



.6J(1), 



a' 



.6J(1), 



a' 



.6J(1), 



a' 



.6J(1), 



a' 



in other words 



yd) ^ 



vS~ 



00c 

\ 5 y 

^E ^ 
Z) 

\ F y 



' ^fc=l 



' -'A;=4 



'^ D^ 
F 

y F y 

^ c^ 

5 

\ A y 



yd) ^ 



' ^fc=5 



= F, 

'^050^ 
A 

\ Cy 

'^ F ^ 
F 

Y F»y 



(88) 



(89) 



(ii). 



(1) 
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f F^ 




f C 0^ 


yi%~ 
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. nS ~ 
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, nS ~ 


BOO 




[o c o) 




\E J 




^ ^y 




^ E 0^ 




^ b\ 




^ D 0^ 


nS ~ 


OOF 


■ >S " 
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. n!i " 
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y L> y 




^A J 




yO E J 



(iii) s 



(1) 



(90) 





f A 0^ 




f D^ 




f^ S ^ 


yt%~ 
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. nS ~ 
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. nS ~ 


coo 




\0 C J 




{F oj 




^ ^y 




^ E 0^ 
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^ F ^ 


n!a ~ 


OOF 


. n'i". ~ 


.4 


. n!i ~ 


D 




^0 D J 




^B J 




^ Fy 



(91) 



The cases s^^' = k/2 and s^^' = 2j are not considered since there they may forbid three different 
real eigenvalues. We will take case (iii) for the following discussion. 
The most general form for a given mass matrix can then be given as 



/ 



\ 



M 



Avi + Ev4 Bv^ + Fvq Dv2 + Cus 

Cf 3 + Dvq Avfj + Ev2 Bvi + Ff 4 

\ Fv2 + Bv^ Cvi + Dv4^ Avs + Evq I 

where vi = {Hi) are the different VEVs of the six Higgs fields present in the model. 



{Mu) 



C'Hl + D^Hl A^H^ + E^Hl B^ H^ + F^ H^ 
\ F^Hl + B^Hl C^Hl + D^Ht A^ Hi + E^ H^^ J 



(92) 



(93) 



(M, 



d)ij 



( A^Hl + F^FJ B^Hl + F^Hl^ D^HJ + C^Hl \ 

C^Hl + D^H^ A^Hl + E^HJ B^HJ+F^HJ 
[f^HJ + B^HI C^Hl+D^HJ A^Hl^E^Hl^ 



(94) 
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f A^m + E^m B^m + F^m D^m + c^m\ 



{Mi)ij ~ 



(95) 



C^'Hl + D^Hl A^Hl + E^HJ B^HJ + F^HJ 

The two light Higgs mass eigenstates which arise by fine-tuning the superpotential shown in 
Eq. m would then be different linear combinations of the six Higgs fields present in the model. For 
a given set of VEVs for the Higgs fields chosen to obtain the desired Yukawa matrices, we can then 
write the light Higgs eigenstates as 

H^,, = Y, ^k=/?;„ (96) 



with <, = (//;,,). 

One may see from the form of the above mass matrices that there is a natural mass hierarchy 
among the up-quarks, down-quarks and leptons which has a geometric interpretation related to the 
position of the different D-branes, as parameterized by the shifts e^. Recall that the Pati-Salam 
gauge symmetry has been broken to the SM by a process which involves splitting the stacks of 
D-branes as shown in Fig. [H The left-handed quarks are localized at the intersections between 
stacks al and b, the right-handed up- type quarks are localized between stacks al and cl, while the 
right-handed down-type quarks are localized between stacks al and c2. Thus, if stack c2 is shifted 
on the torus by an amount €c2 while stack cl is unshifted (ed = 0), then the down- type quark 
masses are naturally suppressed relative to the up- type quarks. Similarly, the left-handed charged 
leptons are localized at the intersection between stack a2 and b, while the right-handed charged 
leptons are localized at the intersection between stacks a2 and c2. Since stack a2 will be shifted 
by some amount ea2, the resulting charged lepton masses will be naturally suppressed relative to 
the down-type quarks. Thus, from purely geometric considerations, we expect 

mu> md> mi. (97) 

It is also clear from the form of the mass matrices that the mass hierarchies between the 
different up-type quarks among themselves, the different down-type quarks among themselves, and 
the different charged leptons among themselves are due primarily to the different Higgs neutral 
component VEVs. Effectively, these values will determined by fine-tuning the superpotential so 
that there are only two Higgs eigenstates as shown in Eq. [96l Ideally, one would like to be able 
to dynamically determine the Higgs eigenstates from first principles rather than fine-tuning the 
superpotential. However, it does not seem possible to do this at least until the issue of moduli 
stablization has been fully addressed. 
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If we define D^ and D^ as the mass diagonal matrices of up and down-type quarks respectively, 
UlMuU^^ = Du, UlMdU^^ = Dd, VcKM = UlUt\ (98) 

where C/* are unitary matrices, then the squared mass matrices are M^MI and MdMjj^. In the 
Standard Model, we can always make quark mass matrices Mu and M^ Hermitian by suitable 
transformation of the right-handed fields. If we take a case that the M^ is very close to the 
diagonal matrix for down-type quark, in other words, Wl and [/^ are very close to the unit matrix 
with very small off-diagonal terms, then 



VcKM ^ U^W"^ ~ [/", 



(99) 



where we have transformed away the right-handed effects and make them the same as the left- 
handed ones. Then the mass matrix of the up-type quarks is 



Mu ~ V^j^jyjDuVcKM- 



(100) 



(ii) An Example 

For different superpartner spectra, we may determine the mass matrices for quarks and leptons 
by running the RGE's up to the unification scale. For example, for tan/3 ~ 50 the CKM matrix at 
the unification scale fi = Mx may be determined to be 



651, 



661] 



VcKM 



^ 0.9754 0.2205 

-0.2203e°-°°3°^ 0.9749 
\ 0.0075e-^9°* -0.0311e^-0°^ 



-0.0026i ^ 
0.0318 
0.9995 



and the diagonal quark mass matrices D^ and D^ may be written 



/o. 



D. 



mt 



\ 



0000139 













0.00404 
1 



\ 



/o. 



Dd = rub 



\ 



00141 







0.0280 

1 



Then the absolute value of M„ turns out to be 

/ 0.000266 0.00109 0.00747 ^ 



\MJ 



mt 



(101) 



(102) 



(103) 



0.00109 0.00481 0.0310 
\ 0.00747 0.0310 0.999 J 

We can then fine-tune the parameters and Higgs VEVs in Eqs. ([93]) . ([M|) . and ([95]) to fit these 
mass matrices. It looks at the first glance that the solution can be easily found, but we should 
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FIG. 10: Ratio of theta functions R^ [er''^'' i93(i7rK(S, e-'^'")]/[?93(0, g-'^'*)] 



keep in mind that the six parameters from the theta function controlled by the D-brane shifts 
and Wilson-line phases are not independent. In doing so, we are essentially constrained by the 
off-diagnonal terms, as can be seen by comparing Eq. [93] with Eq. 11031 For example, the (11) 
element of |M„| is 0.00109, while the corresponding element of the fitted mass matrix is Bvs + Fvq. 
However, the diagonal (33) element of the fitted matrix is given by Av^ + Evq which should be of 
order unity. Thus, to fit the mass matrix we must require 

Bv^ + Fvq 



Avs + Evq 
which effectively means that -j << 1 and ^ << 1. 



0.00109, 



(104) 



In Fig. [T0| we plot the ratio 



(k) 



e-'"^^^^3{mK5,e-'^'^) 



^ 



(5 = 




(105) 



(Ac)=i?3(0,e- 



It can be seen from this plot that the off-diagonal terms will be naturally suppressed relative to 
the diagonal elements. In order to satisfy the conditions shown in Eq. 11041 to fit the up-type quark 
mass matrix, must must allow k to be large enough such that e~'"'^ ~ 0. This in effect means that 
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'd^{i7rK6,e '^'^) = 1, so that 



A 



B 



^7rK(e+l/3)2 ^ _ ,rK(e-l/3)2 



^ ^g-7rK(e+l/6)2^ _g^ g-7rK(e+l/2)2^ ^ 



-7rK(e-l/6)2 



(106) 



On comparing M„ and Eq. (I93p we can obtain a set of parameters which will give the desired 



solution. For example, if we set 



6J(i) 



K = 39.6 and e^^^^ = 0, where e, 



(i)f/ 



,(1) 



e^'^^ 



^6 



e^ ' = 0, and ec 



{l)U 






0, then 

A^ 
B^ 



0.000001, 



0.000266, 
0.236, 



C^ = 0.000001, vl = 0.999, 



D^ = 0.0316, 
E^ = 0.0, 
F^ = 0.0316, 



vl 



0.981, 
0.00481, 



(107) 



vl = 0.0345, 

which may reproduce exactly the mass matrix Eq. 11031 

Similarly, after fixing the Kahler modulus from what is needed to obtain the correct off-diagonal 
terms for the up- type quark mass matrix, we can then obtain a result for M^- For the present case, 
we have ^^ = 1.58 which is obtained from the previous analysis of the soft-terms for a point in the 



mt 



parameter space satisfying all constraints. With the choices e'^^ = 0.061, where ea 



WD _ (1) 



0, 



MD 



.(1) 



0, and e«^ 



.(1) 



0.366 (by Eq. ([85]) ). the desired solution may be obtained with 



A^ = 0.629, 
B^ = 0.0, 
C^ = 0.000098, 
D^ = 0.00158, 
E^ = 0.0, 
F^ = 0.249, 



0.00224, 

0.0, 

1.58, 

0.0, 

0.0445, 

0.0001, 



(108) 



The calculated mass matrix M^ is thus given as 

^ 0.00141 

\Md\ ~ rub 



V 



0.000025 0.000004 ^ 
0.000155 0.0280 0.0 

0.0 0.000000220 1 y 



Dd. 



(109) 



Note that the down-type quark mass matrix and the lepton mass matrix both involve the same 
Higgs fields. Thus, once the parameters needed to fit the down-type mass matrix are fixed, the 
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only freedom in calculating the lepton mass matrix is from the geometric position of each brane. 
To fit the lepton matrix consistent with the down-type quark matrix, we may choose e'^^ = 0.0, 



where €« 



t 9 



0.183, e 



(1)L 



= e« = 0.0, 4^ 


^ = gW = 0.366, then 


A^ = 1.0, 


B^ = 0.000001, 


C^ = 0.000001, 


D^ = 0.0316, 


E^ = 0.0, 


F^ = 0.0316, 



(110) 



so that the fitted mass matrix of leptons is given as 

/ 0.00224 4.74 x 10"^ 4.45 x 10"^ \ 



\M,\ 



rub 



(111) 



/ 



TTT-T 



(112) 



/ 



4.74 X 10"^ 0.0445 2.24 x 10"^ 
\^ 4.45 X 10"^ 2.24 x 10"^ 1.58 

'' 0.00142 3.0 X 10"^ 2.82 x 10"^ \ 

3.0 X 10"^ 0.0282 1.42 x 10~^ 
y 2.82 X 10"^ 1.42 X 10-9 1.0 

We may compare the eigenvalues mT-{0.014, 0.028, 1} for the fitted matrix with the extrapolated 
lepton masses obtained from running the RGEs up to the GUT scale: 

/ 0.000217 ^ 
0.0458 

1 



D, 



m-r 



\ 



(113) 



/ 



Thus, the electron mass comes out to be too big by a factor of six, while the muon mass is 60% 
too small. However, the tau lepton mass does come out correctly. This result can be understood 
by considering that the only difference between the down-type quark and lepton mass matrices is 
an overall exponential factor which results from the additional shift for the leptonic brane, eai- 
Although this result seems to give the wrong answers for the electron and muon masses, it should 
be kept in mind that these are only tree-level results. Indeed, it is of interest that the error in 
the obtained lepton masses seems to increase with decreasing mass. There could indeed be other 
corrections, such as those coming from higher-dimensional operators, which would contribute most 
greatly to the electron and muon masses since they are quite small. 

In short, the above mass matrices can produce the correct quark masses and CKM mixings, and 
the correct r lepton mass at the electroweak scale. The electron mass is about 5 ~ 6.5 times larger 
than the expected value, while the muon mass is about 50 ~ 60% too small. Similar to GUTs, 
we end up with roughly the wrong fermion mass relation nie/m^ = rrid/ms- The correct electron 
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and niuon masses may in principle be generated via high-dimensional operators by introducing 
the vector-like Higgs fields from the ac' sector [67]. Moreover, the suitable neutrino masses and 
mixings can be generated via the seesaw mechanism by choosing suitable Major ana mass matrix 
for the right-handed neutrinos. 

VIII. A COMMENT ON MODULI STABILIZATION 

In the previous section, it has been demonstrated that it possible to obtain Yukawa textures 
which can reproduce the observed fermion mass hierarchies and mixings (extrapolated at the unifi- 
cation scale). Although this is a very interesting result, it is clear from this analysis that the Yukawa 
couplings depend on several parameters which are not determined within the model. Rather, ran- 
dom values for these parameters have been chosen by hand in order to fit the experimental results. 
Thus, it is far from clear that this model in its current form can offer a completely satisfactory 
explanation for the observed fermion mass hierarchies and mixings. 

Given a concrete string model, the low-energy observables such as particle couplings and result- 
ing masses are functions of the open and closed string moduli. In particular, we can see from the 
analysis of the previous section that the Yukawa couplings depend on the position of each stack of 
branes on the tori as well as on the Kahler moduli. Thus, if these moduli can be fixed dynamically, 
then the possible Yukawa mass textures would be tightly constrained. Indeed, once these moduli 
are fixed, the only remaining freedom is in the Higgs sector, namely the specific linear combination 
of the six pairs of Higgs states which must be fine-tuned to produce the two Higgs eigenstates Hu 
and Hd of the MSSM. 

The job of fixing the position of each stack of branes on each torus is equivalent to fixing the 
open-string moduli. D-brane constructions typically have non-chiral open string states present in 
the low-energy spectrum associated with the D-brane position in the internal space and Wilson 
lines. This results in adjoint or additional matter in the symmetric and antisymmetric represen- 
tations unless the open string moduli are completely frozen. These light scalars are not observed 
and the succesful gauge unification in the MSSM would also be spoiled by their presence. While it 
may be possible to find some scenarios where the problems created by these fields are ameliorated, 
it is much simpler to eliminate these fields altogether. One way to do this is to this is to construct 
intersecting D-brane models where the D-branes wrap rigid cycles, which was first explored in [la ] 
and [171] in the context of Type II compactifications on T^/(Z2 x Zg) which is the only known 
toroidal background which posseses such rigid cycles. 
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For T^ /{Tj2 X ^2) t^^ twisted homology contains collapsed 3-cycles. There are 16 fixed points, 
from which arise 16 additional 2-cycles with the topology of P^ = 5^. As a result, there are 32 
collapsed 3-cycles for each twisted sector. A D6-brane wrapping collapsed 3-cycles in each of the 
three twisted sectors will be unable to move away from a particular position on the covering space 
T^, and thus the 3-cycle will be rigid. 

A fractional D-brane wrapping both a bulk cycle as well as the collapsed cycles may be written 
in the form 



n^ 

^^n 



i^^ + i ( E 4.n^J +i ( E on-, J +M E ^t^k^Ka] ■ (114) 

where the Z)6-brane is required to run through the four fixed points for each of the twisted sectors. 
The set of four fixed points may be denoted as S^ for the twisted sector g. The constants e^ ^ •, e^ ■, 
and e^'^j denote the sign of the charge of the fractional brane with respect to the fields which are 
present at the orbifold fixed points. We refer the reader to [l7[ for a detailed discussion of model 
building on this background. 

Let us consider a local super symmetric model consisting of three stacks of D6 branes wrapping 
fractional cycles with the bulk wrapping numbers and intersection numbers shown in Table IXll This 
model is essentially equivalent to the model we have been studying, modulo the fact that the D6 
branes are now wrapping fractional cycles and there are no longer any flat directions for the open- 
string moduli. Thus, the positions of the D-branes are frozen since they are unable to move away 
from the flxed points. Phenomenologically, this is very desirable. First, the existence of the adjoint 
which results from unstablized open-string moduli can destroy the gauge-coupling unification and 
the asymptotic freedom of SU{2>)c- Second, fixing the open-string moduli means that there is 
limited freedom in the values taken by the brane-shift parameters (e^ of the previous section). 
As we have seen, these parameters play a fundamental role in generating the mass hierarchies 
between up-type and down-type quarks, and the leptons. In addition, the open-string moduli must 
be fixed in order to calculate parameters such as /i and tan j3 which play a very important role 
in determining the low-energy superpartner spectra and the Yukawa couplings. Finally, instanton 
induced superpotential couplings may also be very important for addressing issues such as neutrino 
masses, inflation and supersymmetry breaking. The Euclidean D-branes necessary for calculating 
such effects must wrap rigid cycles. 

In order to have such rigid cycles, we must make a choice of discrete torsion which is related to 
the sign of the orientifold planes. Namely, in order for the background to have discrete torsion, there 
must be an odd number of 06^^'+' planes. Thus, the conditions necessary for tadpole cancellation 
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TABLE XI: D6-brane configurations and bulk intersection numbers for a local left-riglit model in Type 
IIA on the T^/(Z2 x ll^) oricntifold, where the D6-branes are wrapping rigid cycles. Although the bulk 
intersection numbers are the same for stacks a\ and a2, it should be understood that the cycles wrapped by 
these stacks go through different fixed points on at least one torus. 





SU{3)c X S[/(2)l X SU{2)r x C/(1)b-l 




iV 


{n\l^)x{n\f)x{n\l'') 


ns 


UA 


b 


6' 


c 


c' 


4, V ij 


f^fe V jk 


4t V kl 


b 
c 


3 

1 
2 
2 


(1,0) X (1,-1) X (1,1) 

(1,0) X (1,-1) X (1,1) 

(1,3) X (1,0) X (1,-1) 

(1,-3) X (0,1) X (1,-1) 






-8 




4 



-3 
-3 


1 
1 


3 
3 




-1 
-1 
2 


-1 
-1 
1 
1 


-1 
-1 
1 
1 


1 
1 
1 
1 



depend directly on the choice of discrete torsion, which will then determine what combination of 
hidden sector branes and/or flux is necessary in order to have a globally consistent model. The 
hidden sector will also be necessary to cancel twisted-tadpoles associated with the orbifold fixed 
points. 

Of course, in the end, our strategy will be to construct the model in the T-dual picure involving 
magnetized fractional D-branes, in the context of flux compactifications. Indeed, it is necessary 
to turn on fluxes in order to stabilize the closed-string moduli. In fact, one such model with an 
equivalent observable sector as the one we have been studying has been constructed as examples of 
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391]. The next step would 



supersymmetric Type IIA Ads and Type IIB Minkowski flux vacua 
then be to combine these two sources of moduli stablization into a single construction with both 
open and closed-string moduli stabilized. Hopefully, we would then be able to uniquely calculate 
both the Yukawa couplings as well as the superpartner spectra. We are presently working on this, 
and hope to report on our progress in the near future. 
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IX. CONCLUSION 

We have analyzed in detail a three-family intersecting D6-brane model where gauge coupling 
unification is achieved at the string scale and where the gauge symmetry can be broken to the 
Standard Model. In the model, it is possible to calculate the supersymmetry breaking soft terms 
and obtain the low energy supersymmetric particle spectra within the reach of the LHC. Finally, 
it is possible to obtain the SM quark masses and CKM mixings and the lepton masses, and the 
neutrino masses and mixings may be generated via the seesaw mechanism. 

Clearly, this model cannot be regarded as being fully realistic until the moduli stablization 
issue has been fully addressed. There are many free parameters which have been fixed in order to 
obtain the desired values for the Yukawa mass matrices and the value of the gauge couplings at the 
unification scale, although it should be kept in mind that these parameters are tightly constrained 
and it is not possible to tune them to just any value. In the case of the Yukawa matrices, the 
free parameters are the Kahler moduli, the brane positions on each torus (open-string moduli) and 
the specific linear combination of states with which have been indentified the two pairs of Higgs 
eigenstates. Although we have chosen specific values for the moduli fields to obtain agreement with 
experiments, it may be possible to uniquely predict these values by introducing the most general 
fluxes. It might also be possible to fix the open string moduli if we require the D-branes wrap rigid 
cycles. However, it seems likely that the Higgs eigenstates would still need to be fine-tuned. 

On the other hand, this does appear to be the first such string-derived model where it is possible 
to give mass to each family of quarks and leptons. Even if we cannot at present uniquely predict 
these values, it must still be regarded as highly significant in that it is possible to come very close 
to getting correct mass matrices and mixings at the unification scale. This suggests that the model 
may be a candidate for a phenomenological description of elementary particle physics in much the 
same way as the MSSM. 

It is also very appealing that the tree-level gauge couplings are unified at the string scale, 
although it is still an open question if the running of the gauge couplings can be maintained all 
the way down to the electroweak scale. The reason for this is that there are chiral exotic states 
present in the spectrum which are bifundamentals under the observable and hidden sector gauge 
groups. We should note that most of these chiral exotic states can be decoupled at the string 
scale and the rest may may be decoupled at an intermediate scale. Even if this were not the 
case, we have found that the hidden sector gauge interactions will become confining at around 
10^ GeV and 10^^ GeV respectively, and so states charged under these groups will not be present 
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in the low-energy spectrum. However, there are exotic states in the spectrum which transform as 
representations of both the hidden and observable sectors which may effect the RGE running of 
the gauge couplings, although we do not expect them to effect the running that much. If we are 
optimistic, then it is possible that these would amount to threshold corrections which might push 
the unification scale up to the string scale. 

If the model does turn out to be a realistic effective description of the observed elementary 
particle physics, then it should be possible to predict the low-energy superpartner mass spectra 
from the model. Besides the D-brane wrapping numbers and closed-string moduli, the superpartner 
spectra depend strongly on the exact way in which supersymmetry is broken. In principle, it should 
be possible to completely specify the exact mechanism, whether through gaugino condensation in 
the hidden sector, flux-induced soft terms, or via instanton induced couplings. In the present 
analysis, we have studied supersymmetry breaking generically via a generic parameterized F-term 
and have shown that it is possible to constrain the phenomenologically allowed parameter space 
by imposing experimental limits on the neutralino relic density and mass limits coming from LEP. 
We have found that the viable parameter space is quite large. Once the experimentally determined 
superpartner mass spectrum begins to take shape it may be possible to find a choice of F-terms 
which will correspond to the observed spectrum. It would be very interesting to explore the collider 
signatures which this model may produce at LHC for the regions of the parameter space which 
satisfy all phenomenological constraints. 

In summary, the model we have studied may produce a realistic phenomenology, once the 
issue of moduli stabilization has been fully addressed. The model represents the first known 
intersecting D-brane model for which mass may be given to each generation of quarks and leptons. 
Furthermore, the supersymmetry breaking soft terms may be studied in the model and may yield 
realistic superpartner spectra. Certainly, the model and the current theoretical tools are not 
presently developed to the point where specific predictions for the known fermion masses and 
superpartner spectrum may be made. However, it is clearly possible for the model to describe the 
known physics of the Standard Model, as well as potentially describing new physics. 
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